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Abstract 

We compute the quantum gravitational contributions to the standard 
model effective potential and analyze their effects on the Higgs vacuum 
stability in the framework of effective field theory. Einstein gravity neces¬ 
sarily implies the existence of higher dimension </> 6 and <f > 8 operators with 
novel couplings /71/2 in the Higgs sector. The beta functions of these 
couplings are established and the impact of the gravity induced contri¬ 
butions on electroweak vacuum stability is studied. We find that the 
true minimum of the standard model effective potential now lies below 
the Planck scale for almost the entire parameter space ( 771/2 (fUt) > 0.01). 
In addition quantum gravity is shown to contribute to the minimal value 
of the standard model NLO effective potential at the percent level. The 
quantum gravity induced contributions yield a metastable vacuum for a 
large fraction of the parameter space in the flowing couplings 771/2- 
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Figure 1: Double logarithmic plot of the NLO effective potential with and without grav¬ 
itational corrections for sample initial values: (1) green/dotted: 771 (mt) = 0.5, ^(mt) = 
0.3 with minimal value V\ min = —4.9 x 10 69 GeV 4 ; (2) red/dashed: 771 (m t ) = 10', 
772 (mt) = 10 9 with V 2 ,min = —4.6 x 10 54 GeV 4 ; (3) purple/dashdotted: rj\ (m t ) = 10 15 , 
772(771 1 ) = 10 2 ° where the minimum has disappeared. (SM) orange: Without gravity, we 
find V SM ,min = -4.0 x 10 110 GeV 4 . 


A central outcome of the recent Higgs boson discovery [1,2| and the absence of new physics 
signals at the LHC is that the standard model (SM) as a quantum field theory stays perturba- 
tively self-consistent all the way up to the Planck scale Mp\ |3,:4]j|. From this perspective the 
conservative scenario of having no beyond-SM-physics up to Mp \—except for gravity—is viable 
and has attracted attention [5 10 Moreover, the measured values for the Higgs pole mass ??Th 
and the top mass m t have an intriguing consequence for the question of stability of the Higgs 
vacuum: The SM lies close to the border of absolute electroweak vacuum stability and metasta¬ 
bility. Vacuum stability is usually studied by determining the renormalization group improved 
effective Higgs potential V((p): If V(4>) develops a negative minimum below the value V ew of the 
electroweak minimum, the SM becomes unstable; if the inverse decay rate for tunneling from 
the false electroweak minimum at </> = 0 ew to the second (true) minimum at d m ; n is larger than 
the lifetime of our universe, then the SM is said to be metastable. The value at the minimum 
kmin = V(0 min ) < Few is very sensitive to the value of nip and m t . State of the art high preci¬ 
sion perturbative calculations |3,4| using the latest experimental data indicate that the SM has 
a negative metastable minimum. However, it lies deep within the Planck regime: The gauge 
invariant value of VAin at next-to-leading order (NLO) precision is (— V^ 0 ) 1 / 4 ~ 10 10 Mi 


1 pi 


11 


While it is fascinating that the SM may be extrapolated to such high energy scales, it is obvious 
that quantum gravitational effects cannot be ignored any longer in these regimes—even for the 
conservative no-new-physics scenario. In consequence, the celebrated statement of metastability 
of the SM based on the above value for V m \ n is spurious. 

This is the motivation for the present analysis. We adopt the conservative viewpoint of having 
no new physics up to M P \ and study the quantum gravity contributions to the SM effective poten- 


x Of course neutrino masses, dark matter and baryogenesis still require some (mild) extensions of the SM. 
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tial. This is done by treating Einstein gravity as an effective quantum field theory 12 . For scales 


below Mpi = 1.22 x 10 19 GeV, the SM coupled to quantum gravity is perturbatively well defined 
albeit of limited predictability due to the necessity of including higher dimensional operators as 
counterterms at every loop order. In fact, the main impact of the gravitational contributions 
to the Higgs effective potential is that non-renormalizability induces higher dimension counter 
terms of the form ^W0 6 and j£r4> 8 into the effective field theory, with a priory undetermined 
couplings rfr and r] 2 . The addition of such higher dimension operators to the tree-level potential 
has been studied in 13,14,11 as a means to parametrize potential new physics effects on top of 
the SM. Here we show that these terms are necessarily present due to the existence of gravity. 
Even if one chooses the new couplings to be absent at low scales, they are turned on at high 
scales by the quantum gravity contributions to the renormalization group equations (RGE). The 
higher dimensional counter terms have a profound effect: For generic positive values of 771 / 2 (rat) 
the true minimum of the SM effective potential is pushed to sub-Planckian scales. It is gener- 
ically metastable or even stable for a large range of initial values of 771/2 at the scale m t . See 
Figure[l]for generic configurations of V{(j>) with and without gravitational contributions. On top, 
quantum gravity effects contribute to the effective potential at NLO with orders of magnitude 
at the percent level. 


Standard model coupled to Einstein gravity. We consider the standard model coupled to 
gravity with a cosmological constant /i0 

£/;!., = 4 + £gb + + v=g (a + <r d.HdM + m 2 \H\ 2 - A |ff| 4 V (1) 

Here k denotes the dimensionful gravitational coupling constant related to Newton’s constant 
k 2 = 327 xG = Mp 2 . As we are interested in the quantum gravitational contributions to the 
Higgs effective potential at the one-loop order, it is sufficient to study the Higgs-gravity sector 
of the above model since the gauge bosons in £ GB and the matter fermions in £ F decouple at 
this leading perturbative order. Moreover, due to the non-renormalizability of gravity we will be 
forced to include higher dimension scalar operators as counterterms 

£ct = V~9 (-Vi k1 \ H \ 6 - V 2 k 4 \H\ 8 ^j , (2) 


carrying their own dimensionless bare couplings 771 and 772 - We seek the one-loop corrections 
to the SM effective potential V{(j>), where we expand the Higgs doublet about the constant real 

background field 0 with H = 4= ( , ^ 1 \ 1 , ) and ipi} E M. The tree-level Higgs potential 

\ (p + + W2 J 


is then given by 


V 


tree 


( 0 ) + ^ 

yrj 2 ^ V 8 ^16 s 


(3) 


2 Note that we could also include an additional non-minimally curvature-coupled scalar term 2 \/—g R p \H\ 2 

into the Lagrangian. However, the freedom of Higgs-held redefinitions and Weyl rescalings may be used to set 

p = 0 in (Til) [15] at the cost of introducing an additional kinetic term. Here we set p = 0 from the start. 
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In the gravitational sector we work in de Donder gauge expanding the metric held as = 
+ k h^ v . This yields the following standard expansions of the Einstein-Hilbert and de Donder 
gauge fixing terms £ g .f.: 

4 V- 9 R + £ g .f. = \Kp P a ^ 8 d 2 h l6 + 0(h 3 ), (4) 

Kj 

2 

KK = 1 + \ + 4 Kg P a0;l5 Ks + 0(h 3 ), 

pafr-yS = ^76 _ YjCx'y ^50 _ ^aS ^ 7 0 ) 

For the quadratic fluctuations of the graviton and Higgs held in our model described by ([Tj) we 
then have 


^»quad 

H “ 


\K0 [P a ^ S (-d 2 + m 2 A ) 


h'yS h 


^2 


a/3 


a 2 + m 2 c 1 <^ 2 - £ \<h [d 2 + - 2 
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with the ehective masses 

m 2 A = ^ ( - m 2 <p 2 + |A0 4 + ^k 2 ?7i0 6 + \Kg 2 (t> 8 + 27l), 

m| = ^ ( - m 2 0 + A (p 3 + |av 2 ?7i0 5 + |k 4 ?7 2 0 7 ), 

mp = —m 2 + 3A0 2 + ^k 2 ?7i(/) 4 + |k 4 ?7 2 (/> 6 , 
m 2 D = -m 2 + A4> 2 + |av 2 ?7i0 4 + |k 4 ?7 2 0 6 . 


(5) 


( 6 ) 


We see that the graviton obtains a small mass ttia generated by the Higgs held. While beyond 
the scope of this letter, it would be important to further analyze the SM coupled to gravity in 
the light of different Higgs mechanisms for the graviton discussed in e.g. 


16,17 


Here we proceed by writing £^ d H = —} 2 VjM !J Vj with the collective quantum held Vj = 
(h^, (pi, ijij). We perform the path integral at 1-loop order to find the gravitational contributions 
F g ^ oop ) = 4W[0] — 4W[(/>] to the ehective 1-loop Higgs potential, where 3 


AV[(p\ = -^/i 4 d J 77^4 ( 9 ln ( p2 + m d) + 3 ln (P 2 - m o) + 1i Kp 2 + "4)G° 2 - m c ) - 4m s]) • 


(2tt 


The relevant dimensionally regularized integral reads 

, rf 4-2e 


h 4 ~ 2£ 


P 1 / 2 2\ 

77 ^- ln(p — m ) = 


(2*; 
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4—2e 


1 

r- + 
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64n 2 e 64n 2 


( m 2 \ 

3 ' 

In — 


L v h J 

~ 2 


(7) 


( 8 ) 


3 We have det M = 64(— d 2 + m A ) 9 [(— d 2 + m A )(d 2 + ml) — 4m^](d 2 + ml) 3 
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with /j, 2 = 47 r/i 2 e The pole in e yields a renormalization of the scalar field couplings m, A, r/i 

and 7/2. Proceeding in the MS scheme gives the following gravitational contribution to the 
renormalized effective potential of the standard model: 


£V[4,} = 


64V mA 


In 


m. 


I 1 


_h +3 TV 

2 / 64V 
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where for conciseness we have definec (3 


C + = x 


m 2 c — m 2 A ± \j ( m 2 C + m A )' 2 — 16 m% 


Av¬ 


is) 


( 10 ) 


The /3-functions. Adding the counterterms necessary to absorb the 1/e-pole terms to the bare 
Lagrangian yields the renormalized Lagrangian. An equivalent statement is that the effective 
potential obeys a renormalization group equation (RGE) 

EEiiA- _ 70 ^eff(0) = 0 , (11) 

where j3\ i = are the /3-functions of the couplings Aand 7 ^ is the anomalous dimension of 
the SM Higgs field. From ([9]) one thus establishes the one-loop /3-functions of the novel couplings 
7/1/2, which take the following form in the {hi 2 m 2 , k 4 A} —> 0 limit: 


/3 (1) 

/3 (1) 


6 V 7 i 1} + tE 

108 A 7/1 - 8A 2 , 

8727L + 

192A 7/2 + 126t ii + |A 2 - 247 7l A 

V_h ?/ 2 - — n 2 - V 

16 t r 2 4 id 4^2 • 


( 12 ) 


The scaling dimension 7 ^ contains the top-Yukawa and electroweak coupling constants yt,gi,g 2 - 
We stress that the A 2 term in [A 1 ^ as well as the A 2 and i]i\ terms in /3^ are quantum gravity 
induced contributions despite the fact that they are not proportional to k, see Figure [2] The 


remaining non-gravitational term in has been reported in 19 , the non-gravitational terms in 


[A^j were considered in 20 , however we disagree with the numerical factor for the Ar/ 2 term given 


there. We hence see that even in the absence of higher dimensional operators ( 7/1 = 7/2 = 0) at an 
initial scale, these terms will be created in the renormalization group flow|^] Including the Higgs 

4 We stress the use of dimensional regularization here, which only sees logarithmic divergences. A (naive) cutoff 
regularization \p\ < Auv of the integral would also induce a AyyTO 2 power divergence in ([8]). However, this term 
does not contribute to the /3-functions in perturbatively renormalized effective field theory. Nevertheless, in the 
framework of the functional renormalization group approach this term contributes, c.f. 18 . 

5 Note that the term under the square root is always positive at <j> = </> min on the parameter space of 771 , 2 - 

6 We assume that contributions to the RG-flow from other higher dimensional operators in the effective field 
theory are consistently neglectable at this order of perturbation theory. It would be important to investigate this 
point in full detail, cf. e.g. 


19 21 
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Figure 2: Diagrams of gravitational processes contributing to the /3-functions of r/i and 
t )2 in (12). Plain lines represent the Higgs field, wiggly lines the gravitons. 


mass and cosmological constant in the analysis leads to quantum gravitational contributions of 
order K 2 m 2 and k 4 T to all the /3-functions including /3\ and /3 m . However, these terms are of 
order 10” 16 or less and absolutely negligible. We therefore set m = A = 0 in the remaining 
analysis. This also puts the electroweak minimum to zero, i.e. h ew = 0. 


Consistent perturbation theory. The renormalization group improved SM effective poten¬ 
tial is traditionally written in the form 


Csm( 0) = A e ff(/r = 0) — , 

with the effective field-dependent coupling constant 

KM) = e 4rw [AM + A$(/<) + A§M 


li=4 


(13) 


(14) 


where T(0) = 7^(/0cilog/n The explicit expressions for the corrections to the effective 
coupling A e fj up to three loops are given in the appendix of [4j. We follow a consistent use of 
perturbation theory along the lines of 122,23jjTT] assuming A ~ h. This is necessary in order for 
the tree-level A0 4 term to receive non-negligible corrections by the one-loop contribution scaling 
as yf ~ h. Perturbation theory in h is applicable, however, it is not the usual loop expansion. 


The expansion of the tree, one- and two-loop contributions of (13) along these lines yields the 
expansion up to next-to-leading order: 

Vsm 0 W = VsmW + dS LO) (0). (15) 

where Vg^ scales as h and 14 m L °^ as ^■ 

Let us now add the gravitational contributions to this picture. In analogy to the SM case 
above we write 

kgrav(0) = heffO = 0)^, (16) 

with the RG-improved effective coupling 

heff(0) =[e 6rw>) ?7i(/i)/=c 2 0 2 + e 8 rw ? 7 2 (/i)K 4 0 4 + 77 $ (0) 


v=4> 


(17) 


Here ,,«(</,) = can be extracted from fol. To obtain a consistent perturbative 

expansion we assume r]in 2 (j) 2 ~ h and r/2^ 4 0 4 ~ h. Expanding V^ rav (0) in h then yields the 
gravitational correction to the SM effective potential up to next-to-leading order K 2 \ 


= CW + b t I LO, (0). 


(18) 
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Figure 3: Plot of Vmin for different initial values r/i (rrit) and r] 2 (mt) at the top mass 
scale (green dots). In the triangular regions the huge ratio yields a breakdown of 

the numerical solution to the RGEs such that no data points are provided. The three 
(yellow) planes fit the data to good accuracy. At the boundaries of the r/i/ 2 -plane, the 
value at the minimum becomes positive (before the minimum disappears). 


Both, Vg^° and have nonzero imaginary parts and in the following we will restrict to the 

real part of the potential, referring to 24 for an interpretation of the imaginary contribution. 


Minimum of the effective potential. The full leading order (LO) potential read^j 

3 ?/ 2 <b 2 

^ L °(0) = V tiee ((j)) - —2/^ 0 4 + ^ 1,92 terms), 


where we suppress the numerically small gauge coupling terms for brevity. The position of the 
true minimum 0 min = /j min is then the second nontrivial solution of A V LO = 0 with flowing 
couplings, guaranteeing gauge invariance of the minimum, c.f. 23,11 . At NLO for the full 
effective potential we then have 


Knn = ^ L °(0min) + + ^av'^ (0min) • (19) 

In Figure [3] we plot the explicit data for K. lin as a function of the initial values and 

V 2 (mt) at the top mass scale. Here we use the initial conditions provided in |4][ll], in particular 
m t = 173.34 GeV and mu = 125.14 GeV, the available higher loop RG equations for the SM and 
our one-loop RGEs for r)i/ 2 ■ In large regions of the parameter space, the value at the minimum 
is well approximated by the three planes depicted in the double logarithmic plot in Figure [3] 
(bounded by the four (red) lines): 


T/ ab ~ 2 

v min — 


V bc ~ 

mm — 


?72tt 4 10 5 ’ 


v cd ~ 

mm — 


r7jK 4 10 5 ' 


7 See 


25 


for a good review on effective potentials in the context of the SM. 
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Figure 4: Portrait of the Higgs vacuum lifetime as a function of 771/2 at the top mass 
scale. The dark meshed (red) region denotes instability; the light meshed (yellow) region 
could not be explored due to numerical issues. The white areas at the boundaries of 
the plot denote the stability region. On the right hand side we magnify the region 

\m/2\ < 0.5. 


Here the lines bounding the planes are parametrized by (a=d) r] 2 = —30?7 2 , (b) r] 2 = 500/7^ for 
771 > 0 , (c) 772 = —771 for rji < 0. Similar values were reported in [lT| for 7/ 2 = 0 and non-flowing 
771 . The minimum disappears approximately on the lines 771 (mt) — 10 14 and r) 2 (mt) — 10 31 
in the 771 / 2 -plane. The value VAin turns positive shortly before, at around — 10 12 and 

f 72 (m t ) ~ 10 30 . Such astronomically high values for 771/2 should be read as measures of scales M 
when new physics arises in the sense of M = Mp\/^fr\[ and M = Mp\/r]^ 1 . In that sense the 
above thresholds reflect the existence of an intriguing stability scale M r-u 10 10 GeV. Finally, for 
positive 77i/ 2 (m t ) < 0.01 the minimum lies beyond the Planck scale. 

Let us briefly discuss the order of magnitude of the NLO gravity contributions to the effective 
potential. The ratio V^^ o ^(0 m i n )/I^ llin evaluated at the minimum ranges between zero and ten 
percent in large regions of the parameter space. At the boundaries, for large initial 771 (m t ) or 
the minimal value Vmin changes sign as indicated above and the relative gravitational 
contribution becomes large. 


Impact on vacuum stability. We have seen that for a huge range of the parameter space of 
?7i/ 2 , the SM coupled to gravity develops a negative minimum Vmin below the electroweak one. 
The lifetime r in units of the age of the universe T v of the false vacuum is usually estimated 
in the SM by 


26 


W = mil V (FWF ex P[TOp] 


In the presence of the dimension six and eight 
counterterms this equation needs to be modified. Numerical studies 27, 28 indicate that r is 


still well approximated by simply replacing |A(/i)| in the exponent of the lifetime expression by 
|A(/i) + 77 e ff(/i)|. In Figure [5] we show the lifetime of the false electroweak vacuum as a function 
of the initial values for i~ji/ 2 at the top mass scale using this lifetime approximation. We see that 
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the SM coupled to gravity generically yields a highly metastable vacuum for a huge part of the 
parameter space in r)y 2 . Notably, the instability only occurs for negative 772 , i.e. 7/2 < —SOrjj. 
Moreover, we find that the ultrashort lifetimes reported in 27 are confined to a very small 


(black) region in 771 / 2 ( 771 *) parameter space once the RG-flow is taken into account. 

Certainly the lifetime formula employed should be taken cautiously as it does not include 
curvature effects [l5,29]. In regions where |Vmin| 1//4 comes close to Mpi the assumption of a 
flat background metric turns inconsistent. In these regions a full analysis expanding the metric 
around curved backgrounds should be employed (cf. e.g. 30 in this context). 

Beyond the instability region derived above the values of the novel couplings 771 and 772 are 
not restricted by present observational data, as so far only single Higgs interactions have been 
probed experimentally (see e.g. 31 for a recent analysis)]^] Therefore also the exponentially large 
values explored here are not excluded and remain perturbative. 

Finally, we note that the reported results on the values of Rnin and the lifetimes are sensitive 
to the values of m t and tuh even at the level of a 2 GeV variation. Mapping this out is left for 
future work. 
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